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We study the two-particle wave function of paired atoms in a Fermi gas with tunable interaction 
strengths controlled by Feshbach resonance. The Cooper pair wave function is examined for its 
bosonic characters, which is quantified by the correction of Bose enhancement factor associated 
with the creation and annihilation composite particle operators. An example is given for a three- 
dimensional uniform gas. Two definitions of Cooper pair wave function are examined. One of which 
is chosen to reflect the off-diagonal long range order (ODLRO). Another one corresponds to a pair 
projection of a BCS state. On the side with negative scattering length, we found that paired atoms 
described by ODLRO are more bosonic than the pair projected definition. It is also found that at 
(/cfci)" 1 > 1, both definitions give similar results, where more than 90% of the atoms occupy the 
corresponding molecular condensates. 
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I. INTRODUCTION 

Recent advancement in the control of Feshbach 
molecules has given rise to many new experimental obser- 
vations in the world of ultracold atomic gas p], 0, H, EL HI • 
At sufficiently low temperatures, fermionic atoms are 
known to form pairs under an attractive interaction. The 
interaction strength can be manipulated by tuning mag- 
netic fields near Feshbach resonance, which is charac- 
terized by a detuning energy between the close channel 
bound state energy and the open channel collision thresh- 
old. A positive detuning leads to a negative scattering 
length, in this regime paired atoms are loosely bound. 
Upon negative detuning, the scattering length becomes 
positive and atoms can form bound molecules, which 
could further condense into a BEC state. Unlike bosonic 
molecules, the statistics of interacting fermionic atoms is 
dictated by Pauli exclusion principle, the ground state is 
thus made up of a large number of modes, even at zero 
temperature. One usually uses a BCS state to approxi- 
mate the ground state at which fermions are paired up 
according to their natural orbits @ . This is very different 
from BEC formed by pure bosons at zero temperature, 
which is well described by a single mode wave function. 
We may however expect, upon a strong enough inter- 
action, paired fermionic atoms become so tightly bound 
that they look just like bosons 0]. In that case, one nat- 
ural question to ask is, how alike are a fermionic pair and 
a boson? In this paper, we address this question by con- 
structing a Cooper pair creation operator and examine 
its bosonic properties across resonance. 

One fundamental feature that distinguishes fermions 
from bosons is the commutation relation between their 
creation and annihilation operators. For bosons, the 
commutator [C, C'~\ is one, while the anticommuta- 
tor {C, C^} is one for fermions. For composite two- 
particles, the corresponding commutator is not exactly 
one [1, [i| [l(| ■ A useful indicator measuring the deviation 
from the bosonic commutation relation is the M-pairs 



normalization factor xm denned by: (0| C M C^ M |0) = 
MIxm, where |0) is the vacuum state. The value of xm 
reflects the correction of Bose enhancement factor, and 
was used to study ground state excitons statistics [1, [§] 
and the connection to quantum entanglement fl(jj |. The 
key quantity was shown to be the ratio Xm+i/xm which 
goes to one for a perfect boson. This ratio will be one 
of our primary indicators of the bosonic characters of 
Cooper pairs. 

However, there has been an ambiguity in defining the 
explicit form of a Cooper pair wave function. Ortiz et al. 
11 1 have given a discussion at length on this matter. In 
12], Yang showed that off-diagonal long range order ex- 
ists in a superconducting state, and is characterized by a 
dominant eigenvector of the two-particle density matrix. 
The eigenvector is sometimes recognized as a Cooper pair 
wave function. On the other hand, the pair projection 
wave function of a BCS state [l3| is also a candidate. 
Both Cooper pair wave functions will be examined in 
this paper. Their bosonic characters are compared and 
we shall discuss their suitability as a bosonic mode in a 
Fermi gas. 

In this paper, we employ the one-channel approach to 
discuss the crossover phenomena at zero temperature [ill . 
EE EH S3 ■ Specifically, a BCS state will be used as our 
ground state wave function 



m=U(u n + v n ai(3i) |0>. 



(1) 



Here a n and f3 n are the annihilation operators of two spin 
components of fermonic atoms, n denotes the quantum 
number of pairing orbit, and u n and v n are amplitudes 
subjected to normalization constraint \u n \ 2 + \v n \ 2 = 1. 
The number of atoms in each spin component is given 



by N = £„ 



In this paper, we will use the solution 



of «„, v n in homogeneous systems. For trapped systems, 
the amplitudes can be determined by methods described 
inRefs. [H,[IH. 
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II. BOSONIC TESTS 

We begin by reviewing some tests on the bosonic char- 
acters of a particle operator. Consider the annihilation 
operator of a composite particle defined by 



C — y \ n a n (} ni 

' * n. 



(2) 



where En \X n \ 2 = 1. The operator C, when sandwiched 
by the BCS state given in (JTJ) , has the following proper- 
ties, 



(° tc ) = |E n A »<««| +EJ A 

([c,ct]) = Ej A «l 2 (l^ 



l 2 KI , 



|5J 2 ) < I- 



(3) 
(4) 



Eq. ([3]) gives the number of composite particles existing 
in the gas. To quantify how 'bosonic' the C molecule is, 
we study the commutator [C, C'~\ . Note however only 
expectation value of the commutator is given in Eq. (|4|). 
not the commutator itself. How close the expectation 
value to unity is a necessary but not sufficient condi- 
tion for C to be bosonic. To actually compare C with a 
pure boson operator, we adopt the views pointed out in 
Refs.fi, S E3- 

It was suggested that the bosonic char- 
acters should be quantified by the normalization factor 
Xm, where 



(0| C M C^ M |0) = M\ XM . 



(5) 



Xm is obviously 1 if C is a perfect boson. It is often more 
convenient to look at the ratio Xm+i/xm, since 



C\M) 



Xm 

XM-l 



M \M — 1) 



cm 



where \£m) is a correction term orthogonal to \M 
and it has the norm given by 



(cmIcm) = 1-M- 



^ +(m _ 1) Xm+x 
Xm-i Xm 



(0) 



(7) 



So we see that the ratio X m+i/xm plays the role of a 
correction of Bose enhancement factor with respect to a 
many body state. It tells us how the gas differs from 
being bosonic, when one more pair of atoms is added 
to or removed from a M-pairs gas. The closer it is to 
one, the less a correction it is. The criterion of a perfect 
boson is xm+i/xm = 1, only then is (€m\£m) zero. In 
this paper we will examine the case with M = N, which 
is the number of atoms of one of the spin components in 
the gas. 

In the case of fermions and if N ^> 1, Combescot et al. 
Q have shown that 



XN+l 
XN 



N 



where zq can be solved from the equation Q 



N 



iA»r 



1 + Zq \X n \" 



(8) 



(9) 



From this last equation the ratio can be solved numeri- 
cally. InRefs. dill, 

the correction factor in ([6|) has been 
studied in exciton systems. Here we will apply Eq. ([9|) to 
atomic Cooper pairs with A„ defined in the next section. 
One of the general features is that the xn+i/xn would 
deviate more from unity when the density of atoms in- 
creases. This is because when the pair density is large, 
atoms within each pair would see the Pauli effect from 
atoms in nearby pairs. 

It is useful to indicate the meaning of Eqs. ([5]) and 
© through a simple model. Let us consider | | 2 = 
(Sfc'<fe c ) _1 f° r W < kc and zero otherwise, where kc is 
the extension of the wave function in momentum space, 
such that the two particle wave function has a spatial 
radius ~ 1/kc- It can then be shown that Eq. ((9]) reads 



XN ~ Z Q V k% 



(10) 



Noting that 1/k^ is the spatial volume of our pair wave 
function, N is the maximum number of pairs in a total 
volume V , the second term in Eq. p0[) is thus the volume 
occupied by all Cooper pairs over the total volume. In 
the BCS limit where a Cooper size is large (kc ~ Fermi 
momentum kp, ~ k F /6ir 2 ), paired atoms are Pauli 
blocked by atoms in between, preventing a bosonization, 
and hence N/zq is nearly zero. While in the BEC limit 
where kc — > oo, each Cooper pair is essentially isolated 
from each other, and this gives N/zq — > 1 [9j. 



III. COOPER PAIR WAVE FUNCTIONS 

We now discuss two choices of A„ in defining the 
Cooper pair wave functions. First, it was shown in 
that the long range correlation (|ri — r^| 



|12|J that the long range correlation {\V\ — v\\ — > oo, 
1 1"2 — 1*2 1 — > oo ) in a paired state is reflected in the eigen 
value decomposition of the two particle density matrix 

G(riy 2 ;r 1)ra ) = (r[)^ (r 2 ) ^ (r 2 ) V« (ri) 

« \^\r m v m n ^i,r 2 )^^r[y 2 ), (11) 

where ip a (r) = J2 n fn (r) «™ and ipp (r) = J2 n fn ( r ) Pn 
are the field operators of the respective species, and {/«} 
is an orthonormal set of natural orbits. The mode func- 
tion d> can be written as 



0l, r 2) 



E„ <£n/n(r 2 )/n(ri) 



y 



(12) 



U'V r , 



in terms of the natural orbits {fn}- 0(ri, r 2 ) is often con- 
sidered as a Cooper pair wave function and 2 J2 n \u n v n \ 
is the number of atoms that condense into Cooper pairs. 
We shall thus construct C with respect to this wave func- 
tion. <^>(ri,r2) is obviously associated with C through 



u n v n 



(13) 
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which gives the explicit expressions 



En \UnV* n \ 2 (\Un\ 2 ~ finf) 

C><? + = L. (15) 

The second term in (C^C) is smaller than unity, while 
the first term is of order N. So as long as the number 
of particles is large, the second term can be dropped. In 
the large N limit, (C'C\ is just the eigenvalue given in 
(fTTjl . We remark that Eq. |13]) was also recognized by 
Leggett |20| as a form of Cooper pair wave function, and 
recently Salasnich et al. have made use of the definition 
to calculate the condensate fraction [lil ]. 

There is however another way of defining a Cooper 
pair based on the single pair projection from a BCS state 
0, HH, EE US ■ By expressing the BCS state as 



i*>= E Eiw 

\ n / j=Q \ n 



uj 2 



C'^ would then be a Cooper pair creation operator. It 
can be shown that C takes the form @, EH, EE E3 : 



C — / , X' n a n (3 n , 

z 'Tl 

V n /U n 



Applying the previous procedures on C", we have 



(17) 
(18) 



(etc") 

<[C",C't]) 



(N) 2 -2(N) + j: n \u:t 

E„ \Vn/u n \ 2 



* Jr. 



(19) 
(20) 



Note that the BCS state given in Eq. (fT6")) is in fact a 
coherent state if C is perfectly bosonic. 



IV. RESULTS IN A UNIFORM GAS 

Before proceeding, let's recap some familiar results in a 
homogeneous BCS gas. The natural orbits are the plane 
wave mode /k (r) = e ik ' r /vi7 and the occupation ampli- 
tudes are given by 



'1± 



2[i 



y/(k 2 - 2^) 2 + 4A 2 



(21) 



where A = Ana (ipp(r)ip a (r)) = —AnaJ^k^kVk is the 
pairing gap, fi is the chemical potential, and a is the 
scattering length. A, a and the atom density p of 




(k a) 



FIG. 1: (Color online) The fraction of condensate particles 
(C*C)/N with two different definitions: (a) A^ ~ u k Vk (solid 
line) (Eq.gi])) and (b)X' k ~ v k /u k (red dashed line) (Eq.(fl9j). 



each species are related by a regularized gap equation 
and a number equation, 



1 

Ana 



P = 



f d 3 k 


1 


J (2tt) 3 


\2y/(k 2 /2 - p) 2 + A 2 


n r 


d 3 k 2 


v~ J 


(2,f N 




(23) 



where the integration can be expressed in terms of el- 
liptic integrals 16]. The Fermi momentum is defined 

1 /3 

as kp — (67r 2 p) , which is the reciprocal of the inter- 
atomic distance. An important dimensionless parameter 
is (kpa) . The BCS limit is indicated by (kpa) 1 <C 
— 1, the BEC limit corresponds to (kpa)~ ^> 1, and the 
crossover occurs at (fcpa) -1 = [H, El[. Some integrals 
used are listed in the Appendix for reference. 

Using Eq. (|2"Tj) for Uk and ik, we evaluate Eq. ([HI 
[T5)l and Eq. (HI [20]). In Fig. Q] we plot the fraction 
of condensation in the g function of the dimen- 

sionless parameter (kpa) -1 . With either choice of Afc, 
the fraction goes to one in the BEC limit (kpa) -1 ^> 1. 
Notice that (C^C) /N is an appreciably higher fraction 
than (C'^C) /N, showing a dominant condensation of 
atoms into the pair wave function defined in (|12|) . 

The expectation value of the commutator ( [C, C'] ) as 
a function of (kpa) -1 is shown in Fig. [2] Again both 
definitions of A^ give unity in the BEC limit, but Eq. 
(fT5| is always closer to one than Eq. (j2"0|) . 

Next we calculate the factor xn+i/xn- By solving 
Eq. JH]) numerically for A& and A' fc , we obtain the ratios 
Xn+i/xn and x'n+i/x'n fr° m Eq. (0). These ratios are 
shown in Fig. [3] as a function of (kpa) -1 . We see that 
Xn+i/xn is closer to one throughout the transition re- 
gion. Together with the tests based on expectation values 
above, A& defined Eq. (fT3]) seems to be a more suitable 
choice for the bosonic description of paired atoms. 

Our calculations indicate an interesting region roughly 
at —2 < (kpa) -1 < 2 where Cooper pairs transit from be- 
ing non-bosonic to bosonic. Note that it does not require 
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FIG. 2: (Color online) Expectation value of the commutator 
(\C, with two different definitions: (a) A^ ~ UkVk (sohd 
line) (Eq. I|15| l) and (b)X' k ~ Vk/uk (red dashed line) Eq. 
JIB). 
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FIG. 3: (Color online) The quantity \n+i/xn with two 
different definitions: (a) Afc ~ fifc^fc (blue solid line) and 
(b)A' fc ~ Vk/uk (red dashed line). 



(fcpa) 1 3> 1 for the emergence of bosonic character. At 
(fc^a) -1 = 1, the fraction of condensation (C'C)/N is 
already 95%, ([C,Ct]) - 0.94 and xw+i/xw - 0.97. 
In particular at the point where the chemical potential 
(j, = ((fci?a) _1 ~ 0.553) [Hj], which is sometimes recog- 
nized as the boundary between BEC and BCS regimes 

[UHi, we have ([ C > Ct ]) = °- 835 ' Xn+i/xn = 0.937. 
The use of definition (118j) gives slightly weaker numbers, 
but a narrower transition region. 



V. CONCLUSION 

To conclude, three indicators were used to quantify 
the bosonic characters of a Cooper pair in an interact- 
ing Fermi gas : (a) (C^C), (b) ([C,C*]), (c) xn+i/xn- 
Two different definitions of a Cooper pair were examined, 
X' n oc v n ju n and A„ oc u* n v n . Our calculations suggest 
that the latter one provides a better description of the 
Cooper pairs as bosonic particles. Moreover, as the frac- 
tion of composite particles goes to one in the BEC limit, 
the gas is basically in its simplest single mode. It appears 
that using either one of the two definitions makes little 
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FIG. 4: (Color online) Momentum distribution (in units of 
fcf) of two different Cooper pair wave functions, Afe ~ iikVk 
(blue line) and A' fc ~ %>k /uk (red dashed line) . In (a) the BCS 
limit, two definitions differ significantly, while in (c) the BEC 
limit, two definitions overlap (see text). 



differences in the strong coupling regime ikpa)^ 1 > 1. 
This is consistent with the results in Ref. [ll|, in which 
the authors addressed the similarity of both definitions. 
As shown in Fig. 4, the difference between (fT3")l and (fTS)) 
on the BCS side is that, the former only takes into ac- 
count of a few momentum states on the Fermi surface, 
whereas the latter one takes the average of all states in- 
side the Fermi sphere. In a weakly interacting gas, only 
atoms lying on the Fermi surface interact effectively, the 
composite particle based on (fl"3|) would thus be more 
bosonic since it takes into account the dominant corre- 
lated states. In the deep BEC limit, the Fermi sphere is 
totally smeared out, and either choice of A„ would weigh 
different momentum states on an almost equal footage, 
resulting in the merge of two different pictures. 

If the system is not in the BEC limit, we have shown 
that the bosonic character of a Cooper pair depends on 
how the pair wave function is defined. Our work here 
is an attempt to identify what definition is more effec- 
tive to reveal a Cooper pair as a boson. From an ex- 
perimental point of view, it is interesting to search for 
observables associated with C or C", so that one can 
probe the the quantum statistics of Cooper pairs directly. 
We also remark that our method can in principle be ex- 
tended to nonuniform systems. However, the coefficients 
u n and v n , which are calculated from the natural orbits 
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[la . Il9j | , do not have closed forms for analytical discus- 
sions. The question how a trapping potential affects the 
bosonic properties of Cooper pairs is an open problem 
for future studies. 



The integrals are expressed in terms of K and E, which 
are respectively the complete elliptic integral of the first 
and second kind. We also need P n , Legendre function of 
the first kind of degree n. 
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Appendix 

We list here some integrals used in this paper. A detail 
analysis can be found in the paper by Marini et al. [l6j 
and the paper by Ortiz et al. [TlJ. In the following list, 
we adopt the following change of variables, 



x Q = jtt/A, x 2 = fc 2 /(2A) 

and introduce shorthand notations 

x 2 - x , E x = + 
1 



§ \l + x /\lxl + l 



-Xq/^J xl + 1. 



(24) 

(25) 
(26) 

(27) 



So we have J k 2 dk = (2A) 3 ^ 2 J x 2 dx. Some integrals that 
appeared in our calculation are listed below 
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/ 
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(2tt) 3 1 

d 3 k 



(2tt) 3 
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2^ Jo 
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(2A) 3/2 
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4tt 2 


UkVk 




(2A) 3 / 2 
8tt 2 


Vk 
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(2A) 3 / 2 


u k 




2tt 2 



,i-~<l.r ( — - 

E x 



x dx 1 



k. 

E x 



El 



v 2 dx- x 



E x + Cx 



(29) 
(30) 
(31) 



1 l \ _ K(k 2 )-2E(k 2 ) 

2 I ~ , 9 , ,n-1/4 



x 2 dx ( — ) = 



E. 



(4 + 1) 



x dx 



IT 







x 



(1 + q)K(n 2 ) - 2qE(n' 
3(o: 2 + I)- 3 / 4 



x 2 dx f 1 

\ E x 

> 771 _ £ 

x 2 dx 

E x + £ x 

2tt (q 2 - 5) P3/2 (q) + 4qPi/ 2 (q) 
35' 



1. 



(33) 
34) 



(X 2 + 1)-V4 



(35) 



Choosing A fe = u k v k / \/J2k' \uk>v k >\ 2 , Eq. © is solved 
with the help of the following integral, 



x 2 dx Try^o + v^o + 1 + 1 / b 

bE 2 + l ~ y/8y/b(b+l) ' 



(36) 



where b is a positive real number. For X' k = 
fe/wic)/VEk' \vk'/u k '\ 2 , Eq. (HJ) is numerically inte- 
grated and solved. 
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